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ABSTRACT. Two cypripedioid orchid genera, Paphiopedilum and Phragmipedium, are listed in Appendix 
I of CITES and are restricted from international trade. Because of their morphological similarity to other 
genera, however, they may be disguised as belonging to one of the other cypripedioids listed along with 
other orchids in Appendix II of CITES. Sequence analysis was performed on the internal transcribed spacer 
region (ITS) of ribosomal DNA of cypripedioid orchids to develop a molecular marker system capable of 
discriminating among rare species in trade. Molecular analyses concentrated on rare cypripedioid orchids 
from the genera Paphiopedilum and Phragmipedium, which are known to be poached from the wild and 
smuggled across international borders disguised as common species. A total of 48 taxa representing two 
genera (Paphiopedilum, N — 43; Phragmipedium, N — 5) have been sequenced and compared for distinc- 
tiveness. Phylogenetic analyses clearly distinguish between these two genera and among other cypripedioid 
genera, with 5—10 fixed nucleotide differences reported between genera. Within a genus, sections of closely 
related taxa are recoverable in phylogenetic analyses, in most cases, with low sequence divergence within 
sections. ITS sequences available in GenBank have been aligned with data generated for this project, 
resulting in a comprehensive sequence library of 151 sequences representing all genera of cypripedioid 
orchids: 70 Paphiopedilum taxa, 16 Phragmipedium taxa, and 14 Cypripedium taxa, as well as represen- 
tatives from Selenipedium and the monotypic genus Mexipedium (Phragmipedium) xerophyticum. Addi- 
tionally, several organelle intron regions have been screened for variation among genera and species. Both 
the chloroplast trnS-M and the mitochondrial NAD1 intron regions, which varied between genera in nu- 
cleotide substitutions and indels, hold promise for increasing ability to distinguish between these orchids. 
The set of DNA markers examined for this project are diagnostic of these genera, appear to be robust, and 
are suitable for rapid assay to avoid unnecessary complication in the legitimate trade of orchids listed in 


CITES Appendix II. 
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INTRODUCTION 


The Orchidaceae is the largest family of flow- 
ering plants, with more than 25,000 species be- 
longing to 800 genera. Import and export of all 
orchids are regulated by the Convention on In- 
ternational Trade in Endangered Species of Wild 
Fauna and Flora (CITES), regardless of specific 
conservation status. The U.S. Department of the 
Interior, Fish and Wildlife Service, Division of 
Scientific Authority is responsible for determin- 
ing the impacts of U.S. trade on orchid species, 
and for developing proposals to amend the 
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CITES appendices. Although ca. 20,000 orchid 
species are not considered at risk, all species are 
regulated because many closely resemble rare 
species of concern. Some rare orchids, particu- 
larly in the cypripedioid (slipper orchids) genera 
Paphiopedilum from Southeast Asia and Phrag- 
mipedium from South America are known to be 
poached from the wild and smuggled across in- 
ternational borders disguised as more common 
species. This smuggling is facilitated by the fact 
that the cypripedioid genera share many mor- 
phological traits through either common descent 
or convergence (Atwood 1984). 

The slipper orchids are a distinct group with 
many primitive features relative to other mem- 
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bers of the family Orchidaceae (Dressler 1993), 
yet they also have been described as advanced 
with specialized structures and on-going speci- 
ation, suggesting that the possession of primitive 
features indicates early divergence during the 
evolution of orchids (Atwood 1984). The five 
genera comprising the subfamily Cypripedioi- 
deae (Cypripedium L., Mexipedium V.A. Albert 
& M.W. Chase, Paphiopedilum Pfitzer, Phrag- 
mipedium Rolfe, and Selenipedium Rchb.f) are 
united by features such as two fertile anthers and 
a saccate lip (pouch) (Averyanov et al. 2003). 
Although several floral differences exist between 
genera, the main vegetative difference is the 
growth pattern. The plicate-leaved genera (Cyp- 
ripedium and Selenipedium) have leaves ar- 
ranged spirally on an elongate stem, where the 
conduplicate-leaved slipper orchids (Mexipe- 
dium, Paphiopedilum, Phragmipedium) have 
leaves arranged in a distichous rosette (Atwood 
1984). Without appreciable vegetative differenc- 
es between species that can be recognized by 
customs officials, it is difficult to accurately 
identify unknown plants that are not producing 
flowers. 

Molecular genetic techniques have been gain- 
ing widespread use in wildlife forensics and can 
potentially be used in U.S. ports to determine 
whether a particular trade specimen is actually 
a Paphiopedilum or Phragmipedium species list- 
ed in CITES Appendix I (commercial trade 
banned), or is simply a species that resembles 
these taxa yet is listed in CITES Appendix II 
and can be legally traded with proof of proper 
cultivation. For genetic techniques to be applied 
to orchid trade, specific genetic characteristics of 
these genera must first be identified, and then 
the appropriate techniques can be disseminated 
to points of entry for the final monitoring of 
specimens being transported. 

A number of DNA markers have been used 
to distinguish between, or determine historical 
relationships between, related orchid species. 
Nuclear ribosomal DNA, mitochondrial DNA, 
and chloroplast DNA, have all been used in past 
studies. One of the most popular gene regions 
used for drawing phylogenetic inferences at the 
generic level and above in plants has been the 
internally transcribed spacer regions (ITS) of nu- 
clear ribosomal cistron (Alvarez & Wendel 
2003). In particular, the ITS region has been use- 
ful in discriminating genera and species of slip- 
per orchids (Cox et al. 1997). Although this gene 
region has contributed to our phylogenetic un- 
derstanding of plants, a review of the molecular 
genetic processes has suggested that the sole use 
of ITS sequences may be misleading with regard 
to phylogenetic inference (Alvarez & Wendel 
2003). Because of a relatively high rate of sub- 
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stitution and maternal inheritance in plants, non- 
coding organelle gene regions have recently 
been used to determine relationships both at 
higher taxonomic levels (e.g., Zhang et al. 2003) 
and among closely related orchid species (e.g., 
Szalanski et al. 2001). Among chloroplast genes, 
the rbcL gene has shed new light on relation- 
ships within the Orchidaceae (Cameron et al. 
1999) and at the genus level among slipper or- 
chids (Albert 1994). Detecting population-level 
polymorphism has been more challenging using 
chloroplast genes, but many are now turning to 
non-coding intron regions that separate coding 
genes throughout the chloroplast genome. Uni- 
versal primers now exist allowing amplification 
of intron regions in a range of higher plants, and 
these regions have been found to be highly poly- 
morphic in mononucleotide repeats (Provan et 
al. 2001). Because of similar modes of inheri- 
tance and utility, the chloroplast genome in 
plants and the mitochondrial genome of animals 
have been considered as natural counterparts for 
use in phylogenetic studies of their respective 
groups (Olmstead & Palmer 1994). The mito- 
chondrial genome of plants, however, is large in 
size and poorly known relative to the animal mi- 
tochondrial genome, with structure and evolu- 
tionary dynamics that are dramatically different 
as well (Palmer et al. 1992). Many aspects of 
the structure of the genome change rapidly, yet 
primary sequence changes very slowly (Palmer 
et al. 1992), and as a result, the genome has not 
been utilized extensively in plant molecular sys- 
tematics. Like the chloroplast genome, however, 
non-coding intron regions are now being utilized 
with universal primers (Demesure et al. 1995), 
and the utility of such regions, especially length- 
change characters, is now being realized in the 
Orchidaceae (Freudenstein & Chase 2001). 

The availability of robust molecular markers 
to discriminate among orchid specimens could 
facilitate commercial trade in species that are not 
an endangered species concern to resource man- 
agers (i.e., those not listed in CITES Appendix 
D. To reach this end, an examination was made 
of variability and identification of discriminating 
characters in multiple DNA regions within and 
between the five genera belonging to the orchid 
subfamily Cypripedioideae. This study describes 
the utility of molecular tools that may aid in the 
taxonomic identification of orchids imported to 
and exported from the United States and to assist 
in effectively enforcing the provisions of CITES. 


MATERIALS AND METHODS 
Sample Collection 


With the help of collaborators, we were able 
to collect tissue from many cypripedioid orchid 
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species belonging to the genera Paphiopedilum 
and Phragmipedium (APPENDIX 1). The speci- 
mens were assigned collection numbers using 
the prefixes: Cp = Cypripedium, Pp = Paphio- 
pedilum, and Pg — Phragmipedium, which ac- 
company the GenBank submissions. These ID 
numbers are useful in distinguishing data unique 
to this study versus the data from GenBank in 
FIGURES 1—4. Two collections of leaf tissue were 
made from the United States Botanic Garden, 
Washington D.C., by the authors and by P. Ford, 
U.S. Fish and Wildlife Service, Arlington, Vir- 
ginia. Additionally, plant tissue collections were 
obtained from P. Simpson, an orchid enthusiast 
and professor at Shepherd College, Shepherds- 
town West Virginia. Several methods of collect- 
ing tissue were attempted, including storing leaf 
punches on ice and flash-freezing leaf-punches 
in liquid nitrogen. The best results for long-term 
sample stability were obtained when leaf tissue 
was crushed into an FTA card (Whatman Inc., 
Newton, Massachusetts) using a pestle; the FTA 
card then was air-dried and stored at room tem- 
perature in individual envelopes. 


Molecular Methods 


Total DNA was isolated from leaf tissue using 
the plant protocol from the PureGene DNA ex- 
traction kit (Gentra Systems Inc., Minneapolis, 
Minnesota), in quantity and quality suitable for 
subsequent laboratory use. The internal tran- 
scribed spacer regions (ITS-1 and ITS-2), sepa- 
rating structural ribosomal DNA genes, were 
amplified from genomic DNA using the poly- 
merase chain reaction (PCR) and universal prim- 
ers (Baldwin 1992). PCR reactions were carried 
out using 1X PCR buffer (10 mM Tris-HCl, pH 
8.3, 20 mM KCl), 2 mM MsCl, 0.2 mM of 
dNTPs, 0.375 uM of each primer, 0.5 mg/ml 
bovine serum albumin, 2.5 U of Taq polymerase, 
and ca. 100 ng of template DNA in a 20 ul re- 
action volume. Thermal cycling used an MJ-Re- 
search (Watertown, Massachusetts) PTC-200 
thermocycler using the following cycle param- 
eters: initial denaturing at 94°C for 2 min.; fol- 
lowed by 35 cycles of 94?C denaturing (1 min.), 
52-56°C annealing (1 min.), 72°C extension (1 
min), ending with a 5 min. extension at 72°C. 
PCR products were purified with Microcon spin 
columns (Millipore Corp., Bedford, Massachu- 
setts). Purified PCR products were sequenced di- 
rectly using amplification primers and conditions 
suggested for Big Dye Cycle Sequencing (Ap- 
plied Biosystems, Foster City, California). Both 
heavy and light strands were sequenced for con- 
 firmation. 

An intron in the chloroplast genome that lies 
between transfer RNA genes for Serine (trnS) 
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and Methionine (trnfM) was amplified using the 
universal primers “‘trnS” and "trnfM" (Deme- 
sure et al. 1995). Since this region is too long 
to obtain overlapping sequence data using only 
the amplification primers (~1700 bp in Quercus 
robur, Demesure et al. 1995), internal sequenc- 
ing primers were developed based upon se- 
quences obtained and aligned with Spiranthes 
romanzoffiana (Forrest et al. 2004, accession 
number AY 363054). The 3’ end of the amplicon 
can be confirmed using the primer Cp2-R1 5'- 
TCCATAATGATGTACCAGAA-3', and the 5' 
end can be confirmed using either Cp2-Fl 5'- 
TTCTGGTACATCATTATGGA-3' or Cp2-F2 
5’-TTGGATTAGTCTTTCTG-3'. 

The mitochondrial nadl intron separating ex- 
ons B and C was amplified and sequenced using 
the primers “nadl exon B" and “nadl exon C" 
from Demesure et al. (1995). Again, the length 
of the amplicon required designing internal am- 
plification primers (—1550 bp in Quercus robur, 
Demesure et al. 1995). Internal sequencing 
primers were designed for the 3' end of the am- 
plicon: nadib-F1 5'-GGGATATACACCA- 
GGGCAAC-3' [20-mer, forward] and nad1b-R1 
5'-GTTGCCCTGGTGTATATCCC-3' [20-mer, 
reverse]; and for the 5’ end: nad1b-F2 5'-CAC- 
CACTTGGGATGGGAAT-3' [19-mer, forward], 
and nadlb-R2 5'-ATTCCCATCCCAAGTGGTG- 
3' [19-mer, reverse]. 

PCR products were sequenced directly using 
the ABI Prism Big Dye Terminator Cycle Se- 
quencing reaction kit utilizing AmpliTaq DNA 
Polymerase, FS (Applied Biosystems, Foster 
City, California). Sequencing reactions were an- 
alyzed by capillary electrophoresis using the 
ABI PRISM-3100 Genetic Analyzer and DNA 
Sequencing Analysis Software (Applied Biosys- 
tems). Forward and reverse sequences were as- 
sembled using Sequencher 4.0 (Gene Codes 
Corporation, Ann Arbor, Michigan). Multiple 
sequence alignment was carried out using Clus- 
tal X ver. 1.4b (Thompson et al. 1994), with gap 
opening/extension penalties of 20/5, respective- 
ly. The resulting alignment was checked by eye 
using MacClade 4.0 (Maddison & Maddison 
2002). 


Molecular Analyses 


Phylogenetic analyses of DNA sequences 
were carried out using PAUP* 4.0b10 (Swofford 
2002). Both distance- and parsimony-based tree 
reconstruction methods were used. Pair-wise Ki- 
mura two-parameter genetic distances (Kimura 
1980) were calculated from the sequence align- 
ment, and the neighbor-joining method was used 
to determine phylogenetic relationships. For par- 
simony analyses, heuristic searches were run 
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TABLE 1. 
are diagnostic of cypripedioid genera. 


Base position 


in alignment 87 102 190 197 235 252 383 543 
Cypripedium C C T T A T G G 
Selenipedium C C T T A GGG 
Mexipedium C AC T G GGG 
Phragmipedium C G T C G GGG 
Paphiopedilum T G C Cl CT G A A 


using unweighted, parsimony-informative char- 
acters with the following settings: starting trees 
for branch swapping were obtained via stepwise 
addition, 10 random addition sequences per run, 
and tree bisection-reconnection (TBR) branch 
swapping on best trees. Alignment gaps were 
treated as missing data and as a 5th base in sep- 
arate runs, and results differed little in topology. 
Bootstrap resampling (Felsenstein 1985) was 
used to assess the support of relationships re- 
covered in phylogenetic analyses. 


RESULTS 


The multiple sequence alignment of 151 ITS 
sequences from cypripedioid orchids both from 
the present study and from GenBank (Cox et al. 
1997) was 909 base pairs (bp) in length. This 
sequence library contains representative samples 
from 57 out of 69 species of Paphiopedilum rec- 
ognized by Cribb (1998), several more recently 
described species (5), and several hybrids, for a 
total of 70 Paphiopedilum species (see APPEN- 
DIX). Included in the data set are 30 taxa that 
have replicate samples, including both Paphio- 
pedilum and Phragmipedium taxa, allowing 
analysis of intraspecific genetic variation for the 
first time in these genera. The nuclear ITS region 
was highly polymorphic, including 699 poly- 
morphic sites and 526 sites that were phyloge- 
netically informative. Sequence divergence 
among genera was generally high (ca. 16-29%), 
with comparisons to Selenipedium being highest 
(28-41%). Within a genus, genetic distances 
were intermediate (2-8%), except between spe- 
cies of Cypripedium, where genetic distances 
between most species were more similar to dis- 
tances between genera (15-2246, data not 
shown). There were 18 diagnostic fixed differ- 
ences (nucleotides) between genera and/or be- 
tween CITES categories (TABLE 1). Also several 
fixed differences were found for sections within 
Paphiopedilum (Parvisepalum, Brachypetalum, 
and Barbata, data not shown). Many fixed dif- 
ferences occurred in alignment gaps, or insertion 
and deletion mutations (indels) between genera. 
For example, a 55 bp deletion was present in 
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Eighteen single nucleotide polymorphisms in the nuclear internal transcribed spacer region (ITS) that 


549 573 616 635 647 676 720 749 751 789 

T/C T CT G G/T T G/T C T A 
C C T G G T A G T G 
T T T A A T G T T G 
T A A A A C G T CT T 
A A T G G T A C/A G A 


most Paphiopedilum relative to the other cypri- 
pedioid orchids, with the exception of the sec- 
tion Parvisepalum species, the most basal sec- 
tion of Paphiopedilum, which retained this 
stretch of bases. 

We analyzed a subset of the aligned sequenc- 
es, (51 from the present study and four from the 
genera Cypripedium, Selenipedium, and Mexi- 
pedium downloaded from GenBank as out- 
groups). Selenipedium chica was used as the 
outgroup for phylogenetic analyses because of 
the high sequence divergence estimates relative 
to other cypripedioid genera and because of the 
basal placement of this genus in a previous phy- 
logenetic analysis of relationships within the 
cypripedioids (Albert 1994, Cox et al. 1997). 
Both neighbor-joining and parsimony analyses 
yielded similar tree topologies (see parsimony 
results, FIGURE 1). Although including gaps as 
characters decreased the number of most parsi- 
monious trees (TABLE 2), the general topology 
remained very similar, so the analysis with high- 
er consistency index (gaps = missing) is shown 
(FIGURE 1). Monophyly of each genus was 
strongly supported in these analyses, as seen in 
the high bootstrap values for these clades 
(100%, FIGURE 1), signifying the ability to dis- 
tinguish between genera using ITS sequence 
data. Generally, taxa belonging to sections with- 
in a genus clustered together and were supported 
by bootstrap analysis (FIGURE 1; sections of Pa- 
phiopedilum following Cribb 1987, 1998). With- 
in sections of Paphiopedilum and Phragmipe- 
dium, distances were <2%, with the exceptions 
of Parvisepalum and Brachypetalum, which 
were slightly higher (2.4—4.6%, FIGURE 1). The 
section Parvisepalum is the most genetically dis- 
tinct Paphiopedilum section, contained the 
greatest genetic diversity within a section (4.6% 
divergence between species, FIGURE 1), and was 
the most basal (1.e., primitive) group of Paphio- 
pedilum orchids based on the ITS sequence data. 
Another interesting result from this analysis is 
the lack of variation within the section Barbata. 
Note the low sequence divergence among spe- 
cies within this section (0.5%, FIGURE 1), which 
is less than the differentiation in the Paphiope- 
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Cypripedium guttatum 

C. yatabeanum 

Mexipedium (Phragmipedium) xerophyticum 
Phragmipedium (Pg) 6-1 boissierianum [I orifolia 
Pg 6-2 boissierianum | (0.56%) 
Pg 18-1 longifolium 

Pg 10-1 /indleyanum i Platypetalum 






Pg 10-2 lindleyanum (0.0%) 
Pg 7-1 caudatum Ñ Phragmipedium 
Pg 17-1 lindenii W (0.42%) 


Paphiopedilum (Pp) 1-1 armeniacum 
Pp 49-1 Vanda M. Pearman Parvisepalum 
Pp 74-1 micranthum (4.64%) 

Pp 80-1 vietnamense 

Pp 33-1 concolor Brachypetalum 
Pp 20-1 niveum (2.38%) 

Pp 70-1 glaucophyllum 

Pp 67-1 victoria-mariae 

Pp 75-1 primulinum 

Pp 75-2 primulinum var.purpurascens f 
Pp 78-1 victoria-regina d 
Pp 53-1 barbigerum 
Pp 10-2 fairrieanum 
Pp 35-1 henryanum 

Pp 56-1 Harrisianum 
Pp 72-1 insigne 

Pp 73-1 insigne 

Pp 79-1 gratrixianum 
Pp 13-1 Airsutissimum 
Pp 13-2 hirsutissimum | 
Pp 13-3 hirsutissimum var. esquirolei | 
Pp 62-1 spicerianum | 


Pp 54-1 glanduliferum 
Pp 63-1 supardii 


g Cochlopetalum 
(0.85%) 






| Paphiopedilum 
(1.79%) 


Pp 38-1 glanduliferum 

Pp 47-1 wilhelminae 

Pp 23-1 randsii 

Pp 50-1 adductum 

Pp 22-1 philippinense 

Pp 76-1 rothschildianum "C. E." 
Pp 61-1 Shillianum 


Pp 37-1 lowii Pardalopetalum 


Coryopedilum 
(1.79%) 


Pp 37-2 lowii var. richardianum : 
Pp 37-3 lowii var. richardianum (0.88%) 
Pp 71-1 hennisianum 
Pp 14-1 javanicum 

Pp 60-1 schoseri 

Pp 65-1 tonsum 

Pp 65-2 tonsum 

Pp 18-1 mastersianum 
Pp 82-1 superbiens 
Pp 81-1 sukhakulii 

Pp 41-1 Ianthe Stage 
Pp 31-1 wardii 

Pp 73-1 bullenianum 
Selenipedium chica 


Barbata 
(0.48%) 





FIGURE 1. Strict consensus of 19,498 most parsimonious trees based upon 296 phylogenetically informative 
sites of nuclear ITS sequence data. Numbers above branches refer to bootstrap support values, with circled 
values for generic and sectional support. Bars to the left of clades indicate sections within genera according to 
Cribb (1987), and numbers under section names are percent sequence divergence within sections. 
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TABLE 2. Variability of gene regions examined and supporting statistics from parsimony analyses. 


Variable 
characters/P.I. 
Gene region analysis characters 
ITS, new data, 925 bp, no gaps 431/296 
ITS, new data, gaps=5" 592/440 
ITS, new + GenBank, no gaps 709/360 
ITS, new + GenBank, gaps=5" 831/419 
CP trnS-trnfM, 1100 bp, no gaps 97/46 
CP trnS-trnfM, gaps=5" 316/101 
Mt nadl, 1559 bp, no gaps 37/18 
Mt nadl, gaps=5" 217/141 
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Taxa Trees Steps 

no. no. no. C.I. 
56 19,498 608 0.7188 
56 84 1108 0.6916 

126 39,790 1192 0.5126 

126 1323 1542 0.3946 
20 2 48 0.9583 
20 3 123 0.9024 
20 10 23 0.7826 
20 2 173 0.8555 


Note: PI. = phylogenetically informative; C.I. = consistency index. 


dilum section Pardalopetalum or the Phragmi- 
pedium section Lorifolia, both of which con- 
tained either individuals from one species or va- 
rieties of the same species. Several species with- 
in section Barbata have been described as a 
“natural hybrid swarm” (P. schoseri, P. javan- 
icum, P. acmodontum; Koopowitz 2000), which 
may help to explain low levels of divergence 
between species. The section Paphiopedilum, 
which received the lowest overall bootstrap sup- 
port among Paphiopedilum sections, appears to 
contain two distinct groups, each of which re- 
ceive higher support: one containing P. hirsutis- 
simum and P. spicerianum, and the other group 
containing the remaining taxa sampled. Al- 
though most sections of Paphiopedilum could be 
distinguished in phylogenetic analysis of ITS 
data, relationships among taxa belonging to the 
sections Coryopedilum and Pardalopetalum 
were an exception. Taxa belonging to the Par- 
dalopetalum group cluster together, but this 
clade is contained within the section Coryope- 
dilum clade. These sections also appeared close- 
ly related in a previous phylogenetic analysis 
(Cox et al. 1997). Few, if any, nucleotide sub- 
stitutions were observed in comparisons of in- 
dividuals belonging to the same taxon, and 
monophyly of these sequences received high 
bootstrap support (FIGURE 1). 

Branch lengths between sections were very 
short (See FIGURE 2C for a comparison of branch 
lengths within Paphiopedilum), signifying few 
differences between sections, yet species within 
sections share polymorphisms (or are genetically 
similar), as they are recovered in both distance- 
based (not shown) and parsimony-based (FIGURE 
1) analyses. Probably as a result of these short 
branches Oe, little genetic variation distinguish- 
ing between sections), relationships between 
most Paphiopedilum sections were not recov- 
ered in analysis of this ITS data set (note the 
collapsed nodes between Paphiopedilum sec- 
tions except Parvisepalum in FIGURE 1). 

A separate phylogenetic analysis combined 


data from this study with that available in 
GenBank (Cox et al. 1997) with the intent of 
maximizing the number of taxa included be- 
longing to the genera Paphiopedilum and Phrag- 
mipedium. 'This analysis included 126 sequences 
representing 15 Phragmipedium and 69 Paphio- 
pedilum taxa, with multiple individuals and/or 
different varieties of some taxa. Mexipedium xe- 
rophyticum (Phragmipedium xerophyticum) was 
used as the outgroup, as this species clearly lies 
outside of the two genera in question based on 
our analysis (FIGURE 1), and the analysis pre- 
sented by Cox et al. (1997). FIGURE 2A shows 
the general structure of the phylogram resulting 
from parsimony analysis using 360 parsimony- 
informative characters (TABLE 2), and note the 
long branch lengths separating genera. All 
Phragmipedium taxa included in this analysis 
were monophyletic (FIGURE 2B, 100% bootstrap 
value). The monophyly of sections Micropetal- 
um and Platypetalum was highly supported, with 
Micropetalum being the most primitive group of 
Phragmipedium, which agrees with results from 
Cox et al. (1997). The section Phragmipedium 
received poorer bootstrap support (67%, FIGURE 
2B), and the Lorifolia taxa clustered together, 
but the monophyly of the group was not sup- 
ported by bootstrap analysis. Where multiple in- 
dividuals of a taxon were included, they gener- 
ally clustered together in this analysis, including 
cases where species formerly considered sepa- 
rate are now synonymized, such as P. sargen- 
tianum Albert & Borge Pett, which is now ac- 
cepted as P. lindleyanum (Rolfe) Rolfe. 
Sections of Paphiopedilum were also recov- 
ered in the combined analysis (FIGURE 2C, 2D). 
Again, Parvisepalum was the most primitive Pa- 
phiopedilum section, and branch lengths leading 
to taxa were long, or in other words, there were 
many substitutions between species. The section 
most closely related to Parvisepalum was Bra- 
chypetalum, which would be expected following 
the designation of these sections in a separate 
subgenus (Brachypetalum) following Cribb 
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_ Mexipedium xerophyticum 515 





Phragmipedium 


Paphiopedilum 


— 5 changes 


FIGURE 2. Results of parsimony analysis on nuclear ITS sequence data for 126 taxa. FIGURE 2A. One of 
39,790 most parsimonious trees based upon 360 phylogenetically informative sites of nuclear ITS sequence data. 
Branch lengths are proportional to amount of change (see scale bar). 


(1987), yet the unity of these sections is not sup- was well supported and nested within the Cor- 
ported by bootstrap analysis. Support for the yopedilum. Included in section Pardalopetalum 
Cochlopetalum section was high, and individu- are two sister species, P. haynaldianum and P. 
als from the same species clustered together. As lowii, which form a well-supported clade within 
in the first analysis, the Pardalopetalum section the section (95% bootstrap value, FIGURE 2C). In 
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FiGURE 2B. 


Subset of phylogeny showing phlogenetic relationships among Phragmipedium species, with 


bootstrap support values above branches and circled values for generic and sectional support. 


section Coryopedilum, controversy exists re- 
garding the status of P. praestans (see Koopow- 
itz 2000). Phaphiopedilum praestans (Rchb.f.) 
Pfitzer has been synonymized with P. glanduli- 
ferum (Blume) Stein (Cribb 1998, Govaerts 
2003), and has been associated with P. wilhel- 
minae, L.O. Williams, all of which appear close- 
ly related in this analysis. The section Paphio- 
pedilum species again formed two clusters (P. 
hirsutissimum and P. spicerianum versus all oth- 
ers), yet in this analysis there is no bootstrap 


support for these clades or for the section. The 
addition of taxa to the Barbata clade in this anal- 
ysis relative to the previous one resulted in a 
substantial decrease in bootstrap support (100% 
in FIGURE 1 versus 65% in FIGURE 2D). Most of 
the intraspecific duplicates in this section did not 
cluster together (e.g., individuals of P. super- 
biens and P. javanicum, FiGURE 2D). Many of 
the Barbata species sequences from GenBank 
had very long branches compared to those in- 
cluded in the first analysis as well (e.g., those in 
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FIGURE 2D.  Subset of phylogeny showing phylogenetic relationships among Paphiopedilum section Barbata 


species, with bootstrap support values as above. 


the clade containing P. venustum and P. apple- 
tonianum, FIGURE 2D). Samples of hybrid origin 
(P. Vanda M. Pearmen, P. Ianthe Stage, P. Shil- 
lianum, P. Harrisianum) clustered with one of 
the parent species (P. delenatii, P. sukhakulii, P. 
rothschildianum, and P. villosum, respectively) 
and not with the other parent. The hybrids did 
not cluster as an intermediate between both par- 
ent species, which were from different Paphio- 
pedilum sections (APPENDIX 1 and FIGURE 2C, 
2D). 


In a preliminary search for novel gene regions 
that can be used to verify results based upon ITS 
sequences, we have screened several Phragmi- 
pedium and Paphiopedilum taxa for variation in 
intron (i.e., non-coding) regions of organelle ge- 
nomes. In the chloroplast genome, we sequenced 
a 1100 bp intron region (£rnS-trnfM) for two 
Cypripedium, five Phragmipedium, and 13 Pa- 
phiopedilum species that, for the latter two gen- 
era, were chosen to represent sections, allowing 
for maximum sampling of potential genetic di- 
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TABLE 3. 
diagnostic of cypripedioid orchid genera. 


Base position 


in alignment T71 94 109 110 115 
Cypripedium A T T A A 
Phragmipedium C C G G C 
Paphiopedilum C T T G C 


Base position 


in alignment 596 605 611 632 644 
Cypripedium A. T C T T 
Phragmipedium G C C T T 
Paphiopedilum G C A G A 


versity. The region contained 97 variable sites, 
of which 46 were phylogenetically informative 
(TABLE 2). Of these informative characters, there 
were a surprising 28 fixed nucleotide differences 
between genera (TABLE 3). These differences in- 
cluded indels, or gaps in the alignment of se- 
quences, which greatly increased the number of 
variable and informative sites to 316 and 101, 
respectively (TABLE 2). Four indel regions (24 
bp) represented fixed differences between these 
genera (TABLE 4). Separate parsimony analyses 
were preformed with gaps treated as missing 
data and as a fifth character. Resulting topolo- 
gies were very similar (see FIGURE 3, from anal- 
ysis with gaps treated as missing). Paphiopedil- 
um and Phragmipedium species were easily dis- 
tinguished using this data set, as can be seen by 
the high bootstrap values and long branch 
lengths between these genera (FIGURE 3). Unlike 
results based upon ITS sequence data (e.g., FIG- 
URE 2B), there was little genetic variation among 
species of Phragmipedium at this chloroplast lo- 
cus (FIGURE 3). Similar to results from the ITS 
analyses, Paphiopedilum micranthum, from sec- 
tion Parvisepalum, was basal to the genus (FIG- 
URE 3). Species from section Brachypetalum 
were intermediate between Paphiopedilum mi- 
cranthum (section Parvisepalum) and other Pa- 
phiopedilum species included in the analysis. 
In a preliminary screening of the mitochon- 
drial nad1 intron, sampling of taxa followed that 
of the chloroplast gene region, where individuals 
were selected to represent sections within 
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Twenty-eight single nucleotide polymorphisms in the chloroplast trnS-trnfM intron region that are 


Phragmipedium and Paphiopedilum, including 
two Cypripedium species, four Phragmipedium, 
and 13 Paphiopedilum species. The multiple se- 
quence alignment also contained Cypripedium 
passerinum and Paphiopedilum delenatii from 
GenBank (Freudenstein & Chase 2001) and was 
1559 bp in length. There were fewer fixed nu- 
cleotide differences between genera at this gene 
region (seven sites, TABLE 5), yet indel regions 
were numerous, including 61 positions that were 
fixed between genera (TABLE 6). A parsimony 
analysis of the nadl data set, with gaps coded 
as a 5th base, included 114 phylogenetically in- 
formative sites and yielded twelve trees (TABLE 
2, FIGURE 4). This nadl gene intron is clearly 
useful for distinguishing between the three gen- 
era represented, with high bootstrap values at 
nodes defining genera (FIGURE 4). Also, there 
was structuring among Paphiopedilum species, 
with the basal species belonging to section Par- 
visepalum and closely related species clustered 
together (e.g., P. insigne, FIGURE 4). The hybrid 
P. Vanda M. Pearman clustered with one of its 
parent species, P. delenatii (FIGURE 4), which 
agrees with results from phylogenetic analysis of 
ITS data (FIGURE 2C). 


DISCUSSION 


The data presented here clearly demonstrates 
the utility of the assayed gene regions for deter- 
mining the genus of unknown slipper orchids. 
Gene regions in the nuclear, chloroplast, and mi- 


TABLE 4. Twenty-four indel positions in the chloroplast trnS-trnfM intron region that are diagnostic of cypri- 


pedioid orchid genera. 


Base position 
in alignment 


Cypripedium 
Phragmipedium 
Paphiopedilum 


262 263 264 265 694 695 696 697 698 699 700 701 798 799 800 801 802 803 804 1002 1003 1004 


T T T T A A T T T A GA C AT A T A T TAG 
A A T T T AG AC AT A T A T TAG 
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FIGURE 3. One of two most parsimonious trees based upon 46 phylogenetically informative sites in an 
analysis of chloroplast trnS-trnfM intron DNA sequences. Bootstrap support values are shown above branches 


with circled values for generic and sectional support. 


tochondrial genomes contained many informa- 
tive substitutions. Data collected from all three 
regions would allow for high accuracy and 
cross-checking of results from any one region, 
although each region had strengths that could be 
drawn upon for certain types of assays. For ex- 
ample, the nuclear ITS region contained the 


highest level of variability and was able to best 
recover relationships among closely related taxa, 
such as sections among species of Paphiopedil- 
um. The best method for assay for the ITS re- 
gion may be DNA sequencing, which would 
yield a high percentage of informative charac- 
ters. The chloroplast trnS-trnfM region had a 
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TABLE 5. Seven single nucleotide polymorphisms in 
the mitochondrial nadl intron region that are di- 
agnostic of genera. 


Base position 


in alignment 398 424 425 426 427 568 6030 


T C T C C C A 
C 
C 


Cypripedium 
Phragmipedium 


TET C TE X C € 
Paphiopedilum A A A A A G 


large number of fixed nucleotide substitutions 
and would likely be assayed effectively using a 
method such as single nucleotide polymorphism 
(SNPs) detection through real-time PCR or per- 
haps a simple fluorometric reaction visible on a 
typical microplate reader. The mitochondrial 
nadl gene intron had the fewest number of nu- 
cleotide substitutions, yet it had a high number 
of indels that appear to be informative between 
slipper orchid genera. In previous studies, indel 
variation has been observed at low levels, such 
as among species of Spiranthes (Chen & Sun 
1998), and at higher taxonomic levels, such as 
within the Orchidaceae (Freudenstein et al. 
1998). To assay these substitutions, DNA se- 
quencing or restriction fragment length poly- 
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morphism (RFLP) analysis will likely be re- 
quired. Once sequence data is gathered, there are 
methods of analysis to consider that treat strings 
of gaps as one evolutionary event, each of which 
can be coded as a separate binary character (see 
Freudenstein & Chase 2001). For example, in 
the nadl data set presented here, the 61 bp with 
fixed indels between taxa likely represent 10 
evolutionary events, or are made up of strings 
of bases 3-18 bp in length (TABLE 6). In this 
way, presence or absence of sequence in one of 
these regions becomes the necessary phyloge- 
netic information. 

The ITS data set with the new data presented 
here combined with that from GenBank contains 
some of the only sequence data where more than 
one individual of a species is represented, and 
several cases point to the importance of includ- . 
ing such replicates. Scattered throughout the ITS 
Paphiopedilum subtree are examples of both 
well-supported clusters of individuals from one 
taxon (e.g., P. lowii, P. glaucophylum, P. hir- 
sutissimum) and cases where individuals of the 
same taxon do not cluster together (e.g., P. 
schoseri, P. javanicum, P. glanduliferum, Fic- 
URE 2C). There are several reasons for why this 
may have occurred. For one, all described Pa- 


TABLE 6. ‘Twenty-four indel positions in the mitochondrial nad1 intron region that are diagnostic of cypripedioid 


orchid genera. 


Base position 


in alignment 433 612 


Cypripedium 
Phragmipedium 
Paphiopedilum 


QQ | 


Base position 


in alignment 682 686 


Cypripedium 
Phragmipedium 
Paphiopedilum 


ol 


Base position 


in alignment 883 


Cypripedium — G C C 

Phragmipedium — G C C 

Paphiopedilum A — — — 
Base position 


in alignment 1249 1250 1251 1265 1266 


Cypripedium 
Phragmipedium 


A A 
Paphiopedilum — — 


ef ar 
ol 


Base position 


in alignment 1422 1423 1424 1426 


Cypripedium 
Phragmipedium 
Paphiopedilum 


A 


T 


614 


687 


C 
C 


1267 
C 


1427 


A 
A 


615 616 680 681 


691 694 


1248 
A 
A. 


1268 1276 1420 1421 
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1273 


>| | 
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FIGURE 4. One of two most parsimonious trees based upon 141 phylogenetically informative sites in an 
analysis of mitochondrial nadl intron DNA sequences. Bootstrap support values are shown above branches with 


circled values for generic and sectional support. 


phiopedilum species may not be good phyloge- 
netic species, instead some designated “‘spe- 
cies" may represent closely related varieties of 
the same species; or in other words, the speci- 
ation process is probably still ongoing in many 
cases, and the species boundaries seen in DNA 
after speciation may not yet have formed. Little 
information is available regarding morphologi- 
cal variation among individuals in the field 


(Koopowitz 2000), and information regarding 
natural intraspecific genetic variation is even 
less known. Morphological intraspecific varia- 
tion may have lead to overzealous *'splitting" of 
some genera into more species than may be nec- 
essary based on genetic differences between 
species (e.g., section Barbata, see Koopowitz 
2000). Given the rarity of plants remaining in 
the wild for many of these species, it becomes 
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complicated to estimate possible gene flow 
among populations. The possibility of other er- 
rors exists also. Some plants may have been 
misidentified because of convergence of a par- 
ticular diagnostic character or coloration. Since 
our methods have been non-invasive (leaf tissue 
was taken from live plants), voucher specimens 
are not available to go back to for verifying spe- 
cies identity. The ideal situation would be to be- 
gin with tissue from herbarium plants with pre- 
served flowers and vegetative material. This, 
however, becomes increasingly difficult when 
working with rare or endangered species. Good 
digital photographs of the plants also may serve 
as a reference in case of future identification 
questions. Sequencing techniques used by dif- 
ferent labs may introduce errors that may be 
enough to cause individuals belonging to the 
same taxon to cluster separately. This would be 
problematic in sections where genetic differenc- 
es between species are few. Identification of 
these types of problems may need to be con- 
ducted on a species by species case. Clearly, ad- 
ditional samples of all taxa need to be sought to 
develop the most extensive baseline data set to 
facilitate the verification of results. These ex- 
amples point to the need for multiple gene re- 
gions to be used along with ITS to distinguish 
between errors in taxonomy, poor resolution of 
molecular data, and to verify results from anal- 
yses based upon a single gene region. Nonethe- 
less, this nuclear gene region clearly allows dif- 
ferentiation between genera and between most 
sections within the genera Paphiopedilum and 
Phragmipedium. 

Additional tissue samples from positively 
identified species belonging to any of the five 
genera of cypripedioid orchids are needed to 
have a comprehensive representation of as many 
taxa as possible. Especially important will be 
representative samples from species in the gen- 
era Cypripedium, Selenipedium, and Mexipe- 
dium for comparisons to be made with Paphio- 
pedilum and Phragmipedium for organelle gene 
regions. We request assistance from all interest- 
ed parties in obtaining tissue from additional 
species of interest. The sequencing protocol we 
have followed to this point appears to allow rap- 
id and accurate identification to genus and often 
to section of the CITES Appendix I genera of 
slipper orchids. 
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APPENDIX. Collection information for slipper orchid samples/sequences included in phylogenetic analyses. 


Genus 


Cypripedium 
Cypripedium 
Cypripedium 
Cypripedium 
Cypripedium 
Cypripedium 
Mexipedium 


Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 


Paphiopedilum 


Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 


Species 


acaule 
guttatum 
parviflorum 
parviflorum 
passerinum 
yatabeanum 
xerophyticum 


acmodontum 
adductum 
adductum 
appletonianum 
argus 
armeniacum 
armeniacum 
barbatum 
barbigerum 
barbigerum 
bellatulum 
bougainvilleanum 
bullenianum 
bullenianum 
callosum 
charlesworthii 
ciliolare 
concolor 


concolor 


dayanum 
delenatii 
delenatii 
dianthum 
druryi 
emersonii 
exul 


Subspecies, variety, 
or cross 


pubescens 
parviflorum 


Authority* 


Sw. 

(Willd.) O.W. Knight 

Salisb. 

Richardson 

Makino 

(Soto Arenas, Salazar 
& Hágsater) V.A. 
Albert & M.W. 
Chase 

M.W. Wood 

Asher 

Asher 

(Gower) Rolfe 

(Rchb.f.) Stein 

S.C. Chen & EY. Liu 

S.C. Chen & EY. Liu 

(Lindl.) Pfitzer 

Tang & ET. Wang 

Tang & ET. Wang 

(Rchb.f.) Stein 

Fowlie 

(Rchb.f.) Pfitzer 

(Rchb.f.) Pfitzer 

(Rchb.f.) Stein 

(Rolfe) Pfitzer 

(Rchb.f.) Stein 

(Lindl. ex. Bateman) 
Pfitzer 

(Lindl. ex. Bateman) 
Pfitzer 

(Lindl.) Stein 

Guillaumin 

Guillaumin 

Tang & FT. Wang 

(Bedd.) Stein 

Koop. & PJ. Cribb 

(Ridl.) Rolfe 


sample 
source** 


Ch, Ch JA LA 4 CN FA 


va Ch Ch Ch Oh PO Ohn Oh Oh OhtaOh Oh Ohh CO 


ON 


NNNDANANAAD 


ID No. 
Ca-01 


Cp-01 
Cp-02 


Pp50-1 


PpO1-1 


Pp53-1 


Pp73-1 


Pp33-1 


CP trnS-fM 
GenB41k 


Mt nadl 
GenBank 


AY 642483 


ITS GenBank 


Z18526 
AY 642767 
AY 642484 
AF3 14830 
Z18527 
Z18515 


Z18446 
AY643459 
278468 
Z18444 
Z78448 
AY643431 
Z18496 
Z18439 
Z18486 
AY643442 
Z18492 
Z18452 
AY643465 
Z18442 
Z18457 
Z18484 
Z18460 
AY643435 


AY642472 


AY642478 


Z18491 


Z18459 
Z18497 
AF314831 
Z18471 
Z18489 
Z18495 
Z18482 


CIC 


VNVAG TS 


SOOT (T'1)9T SUINTOA 


APPENDIX. Continued. 


Genus 


Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 


Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 


Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 


Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 


Paphiopedilum 
Paphiopedilum 


Species 


fairrieanum 
fairrieanum 
fowliei 
glanduliferum 
glanduliferum 
glanduliferum 
glaucophyllum 
glaucophyllum 
godefroyae 
gratrixianum 
gratrixianum 
Harrisianum 


haynaldianum 
hennisianum 
hennisianum 
henryanum 
henryanum 
hirsutissimum 
hirsutissimum 
hirsutissimum 
hirsutissimum 


hookerae 
insigne 
insigne 
insigne 
Ianthe Stage 


javanicum 
javanicum 
kolopackingii 
lawrenceanum 


lowii 
lowii 


Subspecies, variety, 


or Cross 


barbatum X 


villosum 


esquirolei 


rothschildianum X 


sukhakulii 


richardianum 


Authority * 


(Lindl.) Stein 
(Lindl.) Stein 
Birk 

(Blume) Stein 
(Blume) Stein 
(Blume) Stein 
J.J. Sm. 

J.J. Sm. 
(God.-Leb.) Stein 
Rolfe 

Rolfe 


(Rchb.f.) Stein 

(M.W. Wood) Fowlie 

(M.W. Wood) Fowlie 

Braem 

Braem 

(Lindl. ex Hook.) Stein 

(Lindl. ex Hook.) Stein 

(Lindl. ex Hook.) Stein 

(Schltr.) K. Karas. & 
K. Saito 

(Rchb.f.) Stein 

(Wall. ex Lindl.) Pfitzer 

(Wall. ex Lindl.) Pfitzer 

(Wall. ex Lindl.) Pfitzer 


(Reinw. ex Lindl.) 
Pfitzer 

(Reinw. ex Lindl.) 
Pfitzer 

Fowlie 

(Rchb.f.) Pfitzer 

(Lindl.) Stein 

(Asher & Beaman) O. 
Gruss 


Sample 
source ** 


LN Ch Ch Oh ra Ohn Ch Ch tA 


E Ch Li rä Ch th Ch ZG Ch 


— OC 


in 


ON 


Ne HN NO 


ID No. 
Pp10-2 


Pp54-1 
Pp38-1 
Pp70-1 
Pp79-1 


Pp56-1 


Pp71-1 
Pp35-1 


Pp13-1 
Pp13-2 


Pp13-3 


Pp72-1 


Pp72-2 
Pp41-1 


Pp14-1 


Pp37-1 
Pp37-2 


ITS GenBank 


AY 643443 
Z78490 
Z18454 
AY643451 
Z18463 
AY643452 
AY643437 
Z18476 
Z18493 
AY643472 
Z18480 
AY643444 


Z18469 
AY643460 
Z18458 
AY643445 
Z18485 
AY643446 
AY643447 
Z18487 
AY643473 


Z18451 
AY643448 
Z18481 
AY643449 
AY643470 


AY643461 
Z18455 
Z18474 
Z18443 


AY643456 
AY643457 


CP trnS-fM 
GenBank 


AY642470 


AY642471 


AY642473 


AY642480 
AY642475 


Mt nadl 
GenBank 


AY642772 


AY642773 


AY642774 


AY642778 


AY642780 


AY642777 


SONIGHHOO?Id ODOI ANODHS 


EIC 


APPENDIX. Continued. 


Genus 


Paphiopedilum 


Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 


Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 


Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 


Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 


Species 


lowii 


lowii 
malipoense 
mastersianum 
mastersianum 
micranthum 
niveum 
niveum 
papuanum 


parishii 
philippinense 
philippinense 
primulinum 
primulinum 
primulinum 


purpuratum 
randsii 
rothschildianum 
rothschildianum 
sanderianum 
Shillianum 


stonei 
supardii 
supardii 
tonsum 
tonsum 
tonsum 
schoseri 
schoseri 
spicerianum 
superbiens 
superbiens 


Subspecies, variety, 
or Cross 


richardianum 


purpurascens 


‘Charles Edward’ 


(curtisii X 
lawrenceanum) X 
rothschildianum 


braemii 
braemii 


Authority* 


(Asher & Beaman) O. 


Gruss 
(Lindl.) Stein 


S.C. Chen & Z.H. Tsi 


(Rchb.f.) Stein 
(Rchb.f.) Stein 


Tang & ET. Wang 


(Rchb.f.) Stein 
(Rchb.f.) Stein 


(Ridl ex Rendle) L.O. 


Williams 
(Rchb.f.) Stein 
(Rchb.f.) Stein 
(Rchb.f.) Stein 


M.W. Wood & P. Taylor 
M.W. Wood & P. Taylor 
(M.W. Wood) & PJ. 


Cribb 
(Lindl.) Stein 
Fowlie 
(Rchb.f.) Stein 


(Hook.) Stein 

Braem & Lób 
Braem & Lób 
(Rchb.f.) Stein 
(Rchb.f.) Stein 
(Rchb.f.) Stein 


Braem & H. Mohr 
Braem & H. Mohr 


(Rchb.f) Pfitzer 
(Rchb.f.) Stein 
(Rchb.f.) Stein 


sample 
source** 


4 


NO BON bé QN ON 


ZS OS Ota DA 


ANBAR Ch 


Ch ra Lë CNL nh Ohn Cc 


ID No. 
Pp37-3 


Pp18-1 


Pp74-1 
Pp20-1 


Pp21-1 
Pp75-1 


Pp75-2 


Pp23-1 
Pp76-1 


Pp61-1 


Pp63-1 


Pp65-1 
Pp65-2 
Pp60-1 


Pp62-1 
Pp82-1 


ITS GenBank 
AY643474 


Z78472 
Z18498 
AY643466 
Z18449 
AY643432 
AY643436 
Z18494 
Z78450 


Z78470 
AY643475 
Z78466 
AY643438 
Z18479 
AY643439 


Z18440 
AY643458 
Z78465 
AY643453 
Z18473 
AY643471 


Z18467 
AY643454 
Z18475 
Z18456 
AY643463 
AY643464 
AY643462 
Z18453 
AY643450 
AY643467 
Z18441 


CP trnS-fM 
GenBank 


AY642481 
AY642482 


AY642476 


AY642477 


AY 642479 


AY642474 


Mt nadl 
GenBank 


AY 642775 


AY642779 


AY642776 


Vic 


VNVASWTIHS 


SOOT (c'I)Oc 9um[oA 


APPENDIX. Continued. 


Genus 


Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 


Paphiopedilum 
Paphiopedilum 


Paphiopedilum 


Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Paphiopedilum 
Phragmipedium 
Phragmipedium 
Phragmipedium 
Phragmipedium 
Phragmipedium 
Phragmipedium 


Phragmipedium 
Phragmipedium 
Phragmipedium 
Phragmipedium 
Phragmipedium 
Phragmipedium 
Phragmipedium 
Phragmipedium 


Phragmipedium 


boissierianum 


Subspecies, variety, 


Species 


sukhakulii 

sukhakulii 

tigrinum 

urbanianum 

Vanda M. Pearman  bellatulum X 
delenatii 


Or CrOSS 


venustum 
victoria-mariae 


victoria-mariae 


victoria-regina 
victoria-regina 
vietnamense 
villosum 


 wardii 


wardii 

wilhelminae 

besseae 

boissierianum 

boissierianum 

boissierianum 

czerwiakowianum 
caricinum 


caudatum 
caudatum 
caudatum 
caudatum 
exstaminodium 
lindenii 
lindenii 
lindleyanum 


(sargentianum) 


lindleyanum 


(warscewiczianum) 


Authority* 


Schoser & Senghas 
Schoser & Senghas 
Koop. & N. Haseg. 
Fowlie 


(Wall. ex Sims) Pfitzer 


(Sander ex Mast.) 
Rolfe 
(Sander ex Mast.) 
Rolfe 
(Sander) M.W. Wood 
(Sander) M.W. Wood 
O. Gruss & Perner 
(Lindl.) Stein 
Summerh. 
Summerh. 
L.O. Williams 
Dodson & J. Kuhn 
(Rchb.f.) Rolfe 
(Rchb.f.) Rolfe 
(Rchb.f.) Rolfe 
(Rchb.f.) Rolfe 
(Lindl. & Paxton) 
Rolfe 
(Lindl.) Rolfe 
(Lindl.) Rolfe 
(Lindl.) Rolfe 


Castafio, Hagsater & 
E. Aguirre 
(Lindl.) Dressler & 
N.H. Williams 
(Lindl.) Dressler & 
N.H. Williams 
(R.H. Schomb. ex 
Lindl.) Rolfe 
(Lindl.) Rolfe 


Sample 
source** 


NO ra Ch Ch Ch ZS 


Li 


Ch Ch Ch ki bei Ch ba CN bh Ch Ch 


Ch Ch ON LN ka 


UA 


ID No. 
Pp81-1 


Pp49-1 


Pp67-2 
Pp78-1 
Pp80-1 
Pp3 1-1 
Pp47-1 


Pg06-1 
Pg06-2 


Pg07-1 
Pg07-2 


Pg17-1 


Pg10-1 


ITS GenBank 


AY643468 
278462 
Z18488 
Z18445 
AY643434 


Z18447 
Z78477 


AY643440 


Z18478 
AY643441 
AY643433 
Z18483 
AY643469 
Z18461 
AY643455 
Z18513 
AY643424 
AY643425 
Z78502 
Z78503 
Z78510 


AY643429 
Z78501 
Z18500 
Z785 11 
AY 643430 
Z18507 
Z78506 


AY643427 


CP trnS-fM 
GenBank 


AY642467 


AY642465 


AY 642469 


Mt nadl 
GenBank 


AY642781 


AY 642782 


AY 642783 


AY642766 


AY 642770 


SONIGHHOO?Id ODOT ANODHS 


CIC 


APPENDIX. Continued. 








Genus 


Phragmipedium 
Phragmipedium 
Phragmipedium 
Phragmipedium 
Phragmipedium 


Phragmipedium 


Phragmipedium 
Phragmipedium 


Phragmipedium 
Selenipedium 


Species 


lindleyanum 
lindleyanum 
lindleyanum 
lindleyanum 
longifolium 


longifolium 


pearcei 
schlimii 


wallisii 
chica 


Subspecies, variety, 


Or Cross 


(sargentianum) 
(kaieteurum) 


Authority * 


(Lindl.) Rolfe 
(Lindl.) Rolfe 
(Lindl.) Rolfe 


(Warsz ex Rchb.f.) 
Rolfe 

(Warsz ex Rchb.f.) 
Rolfe 

Rauh & Senghas 

(Linden & Rchb.f.) 
Rolfe 

(Rchb.f.) Garay 

Rchb.f. 


Sample 
source** 


5 
5 
6 
6 
5 


ON ON 


ID No. 


Pg10-2 
Pg10-3 


Pg18-1 


CP trnS-fM Mt nadl 
ITS GenBank GenBank GenBank 


AY643428 AY 642468 AY 642768 
Z18505 
Z78503 
AY643426 AY642466 AY 642769 


Z78508 


Z78509 
Z18514 


Z78512 





* Species description authorities and synonyms following Cribb 1998, Govaerts 2003. 


** Sample sources: (1) Shepherd College by P. Simpson; (2) Greenville Co., SC, and (3) Linville, NC, by Lucy Dueck; U.S. Botanic Garden (4) by P. Ford or (5) 
by the authors. Samples from GenBank: (6) nuclear ITS sequences by Cox et al. 1997, with species names used in submission but now synonymized in 
parentheses; (7) mitochondrial NADI by Freudenstein & Chase 2001; (8) ITS sequence for Selenipedium chica by M. Chase pers. comm. 


Note: ID No. 


collection numbers listed in the figures. 


OIC 


VNVATIYS 


SOOT (T’I)9T IUNIO A 


